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Abstract 

We describe a general purpose Mathematical package for com- 
puting Superfield Operator Product Expansions in meromoiphic N — 2 
superconformal field theory. Given the SOPEs for a set of "basic" 
superfields, SOPEs of arbitrarily complicated composites can be com- 
puted automatically. Normal ordered products are always reduced to 
a standard form. It is possible to check the Jacobi identities, and to 
compute Poisson brackets ("classical SOPEs"). We present two ex- 
plicit examples: a construction of the "small" iV = 4 superconformal 
algebra in terms oi N = 2 superfields, and a realisation of the N = 2 
superconformal algebra in terms of chiral and antichiral fermionic su- 
perfields. 
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1 Introduction 



It is well known that Operator Product Expansions (OPEs) form a powerful 
tool in two-dimensional Conformal Field Theory (GET). Once the OPEs for 
the currents are known, it is possible to compute their correlation functions 
algebraically. One also derives differential equations for the other correla- 
tion functions in the theory In gauge theories, one can compute the 
anomaly in the BRST nilpotency, find the spectrum, etc. Although the 
calculations are in principal straightforward. In practice, they can become 
quite cumbersome to do by hand. In the case of a meromorphic GET, these 
calculations can be performed with the OPEdefs package implemented 
in Mathematical'^ |^, an interactive environment for performing symbolic 
computations. 

In the case of super symmetric conformal theories, it is often useful to 
work in a superspace formalism. This drastically reduces the number of 
fields one needs to consider, and also the amount of algebraic manipulations. 
Eor example, a general N = 2 superfield contains 4 components fields (two 
bosonic and two fermionic), therefore one SOPE for the N = 2 superfield 
defines 10 OPEs for the components. Moreover, the checking of one Jacobi 
identity in terms of = 2 superfields is equivalent to 64 identities in terms 
of components. 

Unfortunately, although the formulas to deal with super OPEs are sim- 
ilar to the ones for ordinary OPEs, they are more complicated and calcu- 
lations are more error-prone. This paper describes a Mathematica package 
that enables one to perform computations automatically in the case of iV = 2 
supersymmetry. 

Of course, the simplest super symmetric extension of the conformal field 
theories possesses = 1 supersymmetry. But an A^ = 1 superfield contains 
only one additional fermionic component with respect to the bosonic case 
and so it is possible to use the OPEdefs package. On the other hand, when 
dealing with higher supersymmetry A^ > 3, most of the superfields are con- 
strained, but they can be nicely represented in terms of A^ = 2 ones. In the 
case of A^ = 2, the constraints (chirality or antichirality) can be easily taken 
into account in S0PEN2defs. This is why we decided to limit ourselves for 
the time being to the case of A^ = 2 superfields. 

The package is of a similar nature to OPEdefs. One needs to declare 
the superfields which are used, and give their respective SOPEs. Nonlinear 
algebras can be handled by using a point -splitting definition for composites. 
It is then possible to check the Jacobi identities for the SOPEs. If one 
introduced arbitrary constants in the SOPEs, it is possible in this way to 
find the values for which the algebra is associative, and hence to construct 

(2) 

a super Vl^-algebra. In this way, the superfield version N = 2 Wf' was 
constructed in . Once the basic SOPEs are given, it is in principle possible 
to compute any OPE (or a particular pole) of any composite. Of course, this 

^ Mathematica is a trademark of Wolfram Research Inc. For details, see p|- 
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is restricted by the available memory and CPU-time. The package is able 
to compute SOPEs, but also Poisson brackets ("classical" SOPEs). It has 
been tested by a number of people, already resulting in many publications. 

The package can be extended to perform calculations using dummy in- 
dices. This is still a limited facility. Not all possible simplifications are 
performed when tensors appear which have symmetries when interchanging 
indices. However, for many cases this possibility is quite useful. 

This paper is setup as follows. In a first section we explain our conven- 
tions for N = 2 superspace. Then we derive the formulas for SOPEs. The 
next section consists of a user's guide for the package. We end with two 
explicit examples. 



Notation 

Input for and output from Mathematica is written in typeset font. Input 
lines arc preceded by "In[nj ;=", and corresponding output statements by 
"Out/ijj =", as in Mathematica. 



2 N = 2 superspace. Notations and conventions 

In this section we will introduce our conventions concerning the N = 2 
superfield formalism. 

The N = 2 superspace Z = [z, 9, 9} can be described by one real bosonic 
coordinate z and pair of two conjugate Grassman cordinates 9,9: 

9^ = P = 0, 99 = -99 . (1) 

The integration measure in A'^ = 2 superspace is dZ = dzd9d9 and the 
standard convention for the integration over 9, 9 is assumed. 

To deal with N = 2 superfields '^{Z), it is useful to introduce chiral and 
antichiral spinor derivatives V, V 

which obey the following relations 

V'^ = V^ = 0, {V,v} = -d^. (3) 

Now we can define chiral F and antichiral F superfields which are subjected 
to the following constraints 

VF = VF = 0, (4) 

respectively. 
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It is also instructive to write down the Taylor expansion for a general 
superfield ^{Z) 



$(^i) = E 1 1 + + hiD 



n>0 



1 



•'2 — '712C712: 



^2,2^2 



'J>(^2), (5) 



where 



'12 — [71 — (72, C12 — Pi 



?2, ^12 = ^1 - ^2 + 2 (^1^2 - 02^l) . (6) 



We give the SOPE for the simplest case of = 2 superconformal algebra 
(SCA): 



J(Zi)J(Z2) 



c/4 gl2^12J(^2) 

72 + 72 + 

-^12 ^^12 

euVj{Z2) ei2Vj{Z2) , ^12^125^(^2) 



-Z'12 ■Z'12 

where J{Z) is a general bosonic superfield 



+ 



Z12 



(7) 



3 Formulas for SOPEs 

In principle, the formulas to work with SOPEs follow by expanding every- 
thing in components, applying the standard OPE formulas, and reassem- 
bling again into superfields. However, this approach is impractical. It is 
easier (and conceptually clearer) to construct the formulas by using Taylor 
expansions, contour integration ... in superspace. We will discuss SOPEs 
here and briefly comment on which changes are needed for the Poisson 
bracket case. 

In fact, it will be useful to introduce a notation valid for any number 
of supersymmetries N Q. For a list of numbers [ii], we will write 9^^^ for 
9'^ 9'^ ... and 6l[^l is used for 9'K . . B"-^ . The "complement" N -[i] is defined 
such that = e^-H^W. 

We also introduce projectors V to extract a component of superfield 
or a SOPE. T^g projects on the term without 0's, V\^i^ recovers the term 
proportional to 9^ and so on: 

v.. = [ d^I^l (8) 



We will use the following notation for SOPEs: 

n<=h{A,B) ^^12) 

where h[A, B) is some finite number, and is usually given in terms of the 
conformal dimensions of A and B as h{A, B) = [dA + ds + -^/2]. Note that 
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with this definition (P), depends on the 6\2, but not on zi. The term 

in [AS]o which has no 9i2 dependence is cahed the normal ordered product 
of A and B: 

■.AB:= V^[AB]o. (10) 
With the definition @, we have: 

V^][ABUZ2,ei2) = £^ ^ / de?^ e^i^'^{Z^2T-^A{Zr)B{Z2) , (11) 

where the contour C2 encircles Z2- If we replace the SOPE of A and B by 
the Poisson bracket {A{Zi), B{Z2}pb, we adopt the above formula as the 
definition for [^i?]n with n > in the Poisson brackets ("classical" SOPEs). 

Formulas 

With all these definitions, we can start to derive how we have to compute 
with SOPEs. First, there are some formulas for computing an SOPE of two 
superfields when one of them is a derivative. These follow simply by taking 
derivatives of the SOPE. We give only one example (in = 2): 

[VAB]n{Z2,9\2) = ^^9l2[AB]n-l{Z2,e\2) +Vi[AB]n{Z2,el2) , (12) 

where the subindex on Pi indicates it acts on coordinates Zi. This implies 

Ve[AB]n = V<i[VABU. (13) 
For a term in the regular part of an OPE, we find: 

[ABU = ^[d-AB]o . (14) 
n! 

The SOPE B{Zi)A{Z2) is by analytic continuation equal to A{Z2)B{Zi). 
It remains then to apply a Taylor expansion (^) to the fields in the latter. 
The result is straightforward. We list only the special case for normal or- 
dering, where drastical simplifications occurs because of the projection 
in (0): 

■.BA: = (-l)l^ll^l -.AB: +{-l)\^mY^t}ldiT?^[AB]i (15) 

i>i ^■ 

Applying this formula for A = B a fermionic superfield, one gets an expres- 
sion for '.AA: in terms of the poles in the SOPE A{Zi)A{Z2). 

Associativity of the SOPE means that SOPEs can be computed in any 
order (inside correlation functions). This can be used to derive the Jacobi 
identities. As a first step, we select a particular term in a double SOPE by 
using contour integrals: 

P[,][^P[,][5cy,(Z3) = £^^<3"'^^^iV' 

i, S <3"'''^2V' A{Z,)B{Z2)C{Z,) , (16) 
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where C2 denotes a contour which encircles Z2 once anti-clockwise. We can 
now use a contour deformation argument relating the contour integral in eq. 
( p!6| ) to a contour integral where the integration over Z2 is performed last. 
In fact, the contour deformation works only for the ^-independent part of 
the integral. For the Berezin integral over 0, we just interchange the order 
of integration (taking signs into account). The resulting integral has two 
terms: one where the zi contour is around 2:3, and one where it is around 
Z2- We find: 

P[,][^P[,-][5Cy, = (-l)(l^l+f^l)(^+lil)py][5P[,][^C],],+ 

["9- A v^,^[V^^[ABlCU,^,. (17) 
«>o V / 

This equation has to hold (inside correlators) for consistency of the OPE- 
formalism. It is valid for any integers p, g, i.e. also negative numbers. How- 
ever, in practice we use the Jacobi identities eq. ( |T7| ) for positive p, q as 
equations for the singular part of the SOPEs, while for p ox q zero, they 
define how one should calculate with composites. For example, to compute 
an SOPE with a composite, we insert p = 0, [j\ = into eq. (p^Tj). When 
(7 = 0, [i] = as well, we find: 

■.A :SC::= (-l)l^ll^l :5 : AC:: + : (: A5: -(-l)l^ll^l :5yl:) C : , (18) 



where we used eq. (15). 

For Poisson brackets, the normal ordered product of two operators is 
simply replaced by (noncommutative) multiplication. In particular, this 



means that double contractions (/ < q) in eq. (|17|) for p = 0, [j] = drop 
out. However, the Jacobi identities eq. (p!7|) for positive p,q remain exactly 
the same. 



Improvements 

The rules given in this section up to now are sufficient to compute any 
SOPE, and to reorder any composite into a standard form. For example, to 
compute an SOPE where the first operator is a composite, we could use the 
fact that A{Zi)B{Z2) = B{Z2)A{Zi) and (jl^). However, it is possible to 
construct a rule by using contour integrals to do this in one step. We find: 

n[AB]o C], = \\d'A [BC]i+,], + (-l)l^ll^l l\d'B [AC]i^,], 
/>o'- />o'- 

(19) 

1=1 

where q > 1 and we used eq. (|lj). To have the formula in this simple 
form, we used the convention that any terms with 612 , 612 move to the left 
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of the OPE bracket. For instance, if [BC],n = 9^0, then [A [BC\m\n\ = 
{-l)\''\eu[ADU. 

For composites, we find: 

-.-.AB: C: = :A .BC:: + 

5: i :5'AP0[i3C],: +(-1)1^11^1 E^^'^ P0[AC],: .(20) 

4 User's Guide 

This section is intended as a user's guide to the package S0PEN2defs. Ex- 
plicit examples are given for most operations. S0PEN2defs is based on 
OPEdefs 3.1 and most commands are exactly the same. S0PEN2defs 
requires Mathematica 2.0 or later. 

As S0PEN2defs is implemented as a Mathematica package, it has to 
be loaded before any of its global symbols are used. Loading the package 
a second time will clear all previous definitions of operators and OPEs, as 
well as all stored intermediate results. Assuming that the package is located 
in the Mathematica-path, e.g. in your current directory^, issue: 

In[l] := «S0PEN2defs.m 

After loading S0PEN2defs into Mathematica, help for all the global symbols 
is provided using the standard help-mechanism, e.g. ?OPE. 

Now, you need to declare the operators that will be used. If you want to 
define bosonic (resp. fermionic) operators, use Bosonic (resp. Fermionic). 
For chiral operators, prefix with Chiral, for antichiral ones, use AChiral. 
For instance, to define a bosonic field J and an antichiral bosonic field W, 
the statements are: 

In[2] := Bosonic [J] 

In[3] := AChiralBosonic [W] 

For derivatives we use the following notation: 

d^J ^ J' , DJ^DT[J], PJ^DBT[J], [D, P] J ^ DDB [J] (21) 

The order of the declarations fixes also the ordering of operators used 
by the program: 

J' < DBT[J] < DDB[J] < DT[J] < J < W (22) 

By default, spatial derivatives of an operator are considered "smaller" than 
the operator itself. This can be reversed using the global option NOOrdering 
(see below). 

Finally, the nonregular OPEs between the basic operators have to be 
given by listing the operators that occur at the poles, the first operator in 

■^Use the Mathematica, command SetDirectory for this. 
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the list is the one at the highest non-zero pole, the last operator has to be 
the one at the first order pole. For example, to give J the OPE and W 
of a conformal dimension 1 primary antichiral superfield with C/(l)-charge 
-2: 



In[4] := OPE[J,J] = MakeOPE[{c/4 One + tl2**tbl2** J , 

tbl2**DBT [J] -tl2**DT [J] +tl2**tbl2** J ' }] ; 

In[5] := OPE[J,W] = MakeOPE[{tl2**tbl2**W, 

-2 W - tl2**DT[W]+tl2**tbl2**W'}] ; 

Note the operator One which specifies the unit-operator. The symbols tl2, 
tbl2 denote ^12,^12- In an OPE, they should always be multiplied with an 
operator using **. In particular, one could have tl2**0ne. 

Warning: it is important that the operators occuring as argu- 
ments of OPE in a definition should be given in the order the 



operators are declared (22), otherwise wrong results will be gen- 
erated. 



A normal ordered product : AB : (|10|) is entered in the form NO [A , B] . 
Multiple composites can be entered using only one NO head, e.g. NO [A,B,C] . 
This input is eff'ectively translated into NO [A , NO [B , C] ] . All output is 
normal ordered with the same convention, i.e. from right to left (input can 
be in any order). Also, the operators in composites will always be ordered 



according to the standard order (22). 
In[6] := OPE [J, NO [J, J]] 

Out[6] = « 41 I c tl2**tbl2**0ne/ 4 | |3| | | |2| | c J/2 + 

tbl2**DBT[J] + tl2**DT[J] + 2 tl2**tbl2**N0 [J , J] 
Mill tbl2**DBT[J'] + 2 tbl2**N0[J, DBT [ J] ] + 
tl2**DT[J'] - 2 tl2**N0[J, DT[J]] + 
2 tl2**tbl2**N0[J' , J] » 

This also shows what the output for an OPE looks like. For an OPE with 
n poles, the output is « n| I ... I |n-l I I I 1 1 1 I . . . ». 

Warning: when computing OPEs with composites, or when re- 
ordering composites, S0PEN2defs remembers by default some 
intermediate results. Thus, it is dangerous to change the defi- 
nition of the basic OPEs after some calculations have been per- 
formed. For example, consider a constant a in an OPE. If calcu- 
lations are performed after assigning a value to a, the interme- 
diate results are stored with this value. Changing a afterwards 
will give wrong results. 

The other globally defined functions available from the package are: 

• OPEOperator [operator_ , parity_] provides a more general way to 
declare an operator than Bosonic and Fermionic. The second argu- 
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ment is the parity of the operator such that (— l)P^^*y is +1 for a 
boson, and —1 for a fermion. It can be a symbohc constant. This can 
be used to declare a 6c-system of unspecified parity. In such cases, the 
operator can contain a named pattern: 

In[7] := OPEOperator [J [i_] , parity [i] ] 
If one wants to declare more operators, one can group each operator 
and its parity in a list: 

In[8] := OPEOperator [{b[i_] .parity [i] }, {c [i_] ,parity[i]}] 
Sec also SetOPEOptions [ParityMethod, _] . 

• ChiralQ [operator_] tests if operator is chiral. 
One can use ChiralQ [A] = True. 

• AChiralQ [operator.] is the same for antichiral operators. 

• OPEPole [n_] [ope_] gets a single pole of an OPE: 

In[9] := OPEPole [2] [Out [6] ] 

Out[9j = c J/2 + tbl2**DBT[J] + tl2**DT[J] + 
2 tl2**tbl2**N0[J, J] 
OPEPole [n_] [A_,B_] can also be used to compute only one pole term 
of an OPE: 

In[10] := y. - OPEPole [2] [J, J] //Expand 

Out[10] =0 

OPEPole can also give terms in the regular part of the OPE: 
In[ll] := OPEPole [-1] [J, J] 
Outflij = NO [J' , J] + tbl2**N0[DBT[J'] , J] + 

tl2**N0[DT[J'] ,J] - tl2**tbl2**N0[DDB[J'] , J]/2 

• MaxPole [ope_] gives the order of the highest pole in the OPE. 

• OPEParity [A_] returns an even (odd) integer if A is bosonic (fcrmionic). 

• OPESimplif y [ope_ , Function -> function_] "collects" all terms in 
ope with the same operator and applies function on the coefficients. 
The default setting for the option Function is Expand, but this can 
be changed using SetOptions. 

OPESimplif y [pole_ , Function -> function.] does the same sim- 
plifications on sums of operators. 

The alternative syntax OPESimplif y [ope_ , function.] is also al- 
lowed. 

• OPEMap [function. , ope_] maps function to all poles of ope. 

• OPEMap At [function., ope_, position,] maps function to a part 
of ope specified by position (see MapAt). For example, 

OPEMapAt [Expand , ope , { 1 }] 
maps Expand on the highest order pole of ope. 
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• GetCoef f icients [expr_] returns a list of all coefficients of operators 
in expr which can be OPEs or poles. Applied on a list, it maps on the 
elements of the list. 

• GetOperators [expr_] returns a list of all operators in expr which 
can be OPEs or poles. Applied on a list, it maps on the elements of 
the list. 

• OPEJacobi [opl_ , op2_ , op3_] computes the Jacobi-identities (0) for 
the singular part of the OPEs of the three arguments. In general, all 
different orderings should be tried to ensure associativity. 
OPEJacobi returns a list of which all should be zero up to null fields 
to be associative. 

OPEJacobi accepts an option Function -> function which it passes 
to OPESimplify at intermediate steps. 

• Delta [i_,j_] is the Kronecker delta symbol 5ij. 

• Epsilon[i_, j ] is the antisymmetric symbol (with any number of 

arguments. 

• N20PEToComponents [J_] gives the = components of a superfield 
J. For J = A + eB + eC + eOD the components are {A, B, C, D}. In 
S0PEN2defs, the components are represented as 

{J,VJ,VJ,-^[V,V]J}, 

where a projection on the 9 independent part is understood. 
N20PEToComponents [ope_ , Jl_ , J2_] computes the 16 OPEs of the 
components of Ji and J2- It is a double list where the (m, n)-th ele- 
ment is the A = OPE of the m-th component of Ji with the n-th 
component of J2. 

• ClearOPESavedValues [] clears all stored intermediate results, but not 
the definition of the operators and their OPEs. To clear everything, 
reload the package. 

• TeXForm [ope_] gives output for an OPE. The arguments are 
always Zi , Z2 . 

• TeXFormTD can be assigned a string which will be used by tt TeXForm 
for the output of DT, DBT, DDB. 

• OPESave [f ilenanie_] (with filename a string between double quotes) 
saves the intermediate results that OPEdefs remembers to file (see the 
option OPESaving below). 

• SetOPEOptions is a function to set the global options of the package. 
The current options are: 
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— SetOPEOptions [NOOrdering, n_] : if n is negative, order higher 
spatial derivatives to the left (default), if n is positive, order them 
to the right. 

— SetOPEOptions [ParityMethod, 0|1] : makes it possible to use 
operators of an unspecified parity. When the second argument 
is (default), all operators have to be declared to be bosonic or 
fermionic. When the argument is 1, OPEOperator can be used 
with a symbolic parity. Note that in this case, powers of —1 are 
used to compute signs, which is slightly slower than the boolean 
function which is used by the first method. 

This option is not normally needed as the use of OPEOperator 
with a non-integer second argument sets this option automati- 
cally. 

— SetOPEOptions [OPESaving, boolean_] : if boolean evaluates 
to True (default), OPEdefs stores the intermediate results when 
computing OPEs of composites and when reordering composites. 
This option is useful if Mathematica runs short of memory in a 
large calculation, or when computing with dummy indices (see 
section |5.2| ). 

— SetOPEOptions [OPEMethod, method_] : with method set to Quan- 
tumOPEs enables normal OPE computations (default setting), 
while ClassicalOPEs enables Poisson bracket computations. Us- 
ing this option implicitly calls the function ClearOPESavedValues [] . 

— SetOPEOptions [EnableDummies , True] : loads the Dummies' 

package and sets it up. 

5 Examples 

5.1 The = 4 superconformal algebra 

As a first example of using the SOPEdefs package, let us construct the 
OPEs of TV = 4 SU{2) SCA § in terms of iV = 2 superfields. The iV = 4 
SU{2) SCA contains the following component currents: the stress-tensor 
T{z), affine su{2) currents J^{z),i = 1,2,3 and two doublets of fermionic 
spin 3/2 currents Ga, G^, a, (3 = 1,2 which form respectively the fundamen- 
tal and conjugated representations of su{2). Due to the spin structure, the 
simplest way to put these component currents into N = 2 supermultiplets 
is to start with a general spin 1 = 2 superfield J{Z) which forms the 
N = 2 SCA (spin contents is {1, 3/2, 3/2, 2} together with two bosonic spin 
1 chiral-antichiral superfields W{Z),W{Z), VW = VW = (spin content 
{1,3/2,1,3/2}). 

First of all we need to define all superfields we are dealing with: 
In[l] := «S0PEN2defs.m 
In[2] := Bosonic [J] 
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In[3] := AChiralBosonic [W] 
In[4] := ChiralBosonic [WB] 
As the next step we define the OPE between all the super currents: 
In[5] := OPE[J,J] = MakeOPE[{c/4 One + tl2**tbl2** J , 

tbl2**DBT [J] -tl2**DT [J] +tl2**tbl2** J ' }] ; 
In[6] := OPE[J,W] = MakeOPE[{tl2**tbl2**W, 

-2 W - tl2**DT[W]+tl2**tbl2**W'}] ; 
In[7] := OPE[J,WB] = MakeOPE[{tl2**tbl2**WB, 

2 WB+tbl2**DBT[WB]+tl2**tbl2**WB'}] ; 
In[8] := OPE[W,WB] = MakeOPE[{al tl2**tbl2**0ne, 

a2 One + tl2**tbl2**J , a3 J + a4 tl2**DT [J] }] ; 
Let us note that due to the first OPE, J{Z) obeys the N = 2 SCA, while 
the second and third OPEs define W{Z) and W{Z) as antichiral-chiral su- 
perfields primary spin 1 with respect to J{Z). As concerning the last OPE 
we introduce the arbitrary coefficients al — a4 for all the terms because we 
want to determine the exact form of this OPE. The possible terms and their 
structure in this OPE are completely fixed by the antichiral-chiral structure 
of the W, W supercurrents (the coefficient before term 6*12 * *^i2 * *J can be 
chosen to be 1 due to the scaling invariance W — > aW). 

Now we will fix all coefficients from the Jacobi identities 0: 
In[9] := OPEJacobi [J,W,WB, Function -> Factor] //Union 
Out[9] = {0, ((2-a3)*J)/4, ((2-a3)*J)/2, ((-2+a3)*J)/2, 
(a3-a4)*J, ( (al-a2) *One) /2 , (al-a2)*0ne, 
(-al+a2)*0ne, ((4*al+c)*0ne)/4, 
( (8*a2+a3*c) *One) /4 , ( (8*a2+a4*c) *One) /4 , 
( (-2+a4) *DT [J] ) /2 , (-a3+a4) *DT [J] } 
We can solve the system of equation inside Mathematica, e.g. 
In[10] := soil = Solve [GetCoefficients[7,]==0] 
Out[10]={{al -> -c/4, a2 -> -c/4, a3 -> 2, a2 -> 2}} 
Now we can check that all other Jacobi identities are satisfied: 

In[ll] := OPESimplify [OPEJacobi [W,W,WB]/. First [soli]] //Union 
Out[ll] = {0} 

By this we finished the construction of the = 4 SU (2) SCA in terms of 
N = 2 superfields. 

5.2 A Miura-like realisation of = 2 

As a second example we will demonstrate how to use the S0PEN2defs pack- 
age to find a realisation of the N = 2 SCA in terms of M pairs of chiral- 
antichiral fermionic superfields (Miura type realisation). We also show here 
our technique to deal with summing over repeated indices, which is imple- 
mented in a separate package Dummies. This package can also be used 
together with OPEdefs, or even independently. We only explain some basic 

^ This statement takes 5 seconds on a Pentium 90 MHz running Mathematica 2.2 for 
Windows. 
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commands of Dummies here. 

The first step is loading the package. 

In[l] := «S0PEN2defs.m 
Enable working with dummy indices. 

In[2] := SetOPEDptions [EnableDiimmies , True] 
We now define which indices we are going to use. Ah summation indices 
win have the form i [1] , i [2] , etc. That is, they ah have Head i. 

In[3] := Def ineDummy [i] 
For convenience, we define an abbreviation for the summation range. 

In[4] := dimension [i] = M; 
Now we define the superfields we are going to play with. In the case at 
hand we can define M pairs of chiral-antichiral fermionic superfields Fi , Fj 
through two Mathematica statements: 

In[5] := ChiralFermionic [F [_] ] 

In[6] := AChiralFermionic[FB[_]] 
Using indices the SOPEs between superfields Fi and F^ can be defined 
mostly as they are defined in the textbooks: 

In[7] := OPE [F [i_] , FB [ j _] ] : = 

Delta[i,j] MakeOPE [{-1/2 tl2**tbl2**Dne , One}] 
Now we will introduce a composite superfield J{Z) which should satisfy the 
N = 2 SCA eq. (0): 

J = aiJ2- FiF" ■ +a2VFi + a^VF^ (23) 

i 

where besides the standard term -FjF* we also introduced the N = 2 analog 
of a Feigin-Fuchs term. In Mathematica this becomes: 
In[8] := J := NewDummies [ 

al NO[F[i[l]] ,FB[i[l]]] + 

a2 DBT[F[1]] + a3 DT[FB[1]]] 
Note the NewDummies statement and the assignment using ':='. Together, 
they make sure that everytime J is used, the summation indices will have a 
new number. 

After these definitions we are ready to do real calculations. For example. 



to check for which values of parameters ai , a2 , 03 the superfield J (^31) spans 
N = 2 SCA, let us calculate the SOPE J with J and subtract the known 
result 0: 

In[9] := OPESimplify [DummySimplifyC 
OPE [J, J] - 

MakeOPE [{c/4 One + tl2**tbl2** J , 

tbl2**DBT [J] -tl2**DT [J] +tl2**tbl2** J ' }] 
] , Function -> Factor] 
Outl9] = « 21 I One (-c + 4 M al-2 - 8 a2 a3)/4 + 

(-1+al) a2 tl2 ** tbl2 ** DBT[F[1]] + . . . | | 1| | 
(-1+al) al tbl2 ** NO [DBT [F [i [1] ] ] , FB[i[l]]] + 
...» 
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where DummySimplif y renumbers the dummy indices to get a unique form^ 
For brevity, we omitted a number of terms (all proportional to oi — 1) in 
the output and replaced them by an ellipsis. As can be immediately seen, 
the coefficient ai must be equal to 1, while 02,03 are arbitrary. Let us note 
that our result Out[9] also gives the expression for the central charge of the 
N = 2 SCA: 

c = 4M - 80203 (24) 

Thus we reproduced the standard Miura realisation for the N = 2 SCA. 

We hope that these simple examples will help users to use the S0PEN2defs 
package in more complicated calculations. 

6 How to get it, and the future 

If you are interested in S0PEN2defs, you can get it by Email from the 
authors. You'll also find it at http://euclid.tp.ph.ic.ac.uk/~krthie/ 
or via anonymous ftp at euclid.tp.ph.ic.ac.uk. Please put a reference 
to this paper in your paper when you use it. Questions, remarks and im- 
provements are welcome. The package wil be extended to other numbers of 
supersymmetries A^. 
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